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Abstract 

We calculate the polarisation of ifi's and T's produced by the fragmentation of heavy- 
quarks. We find that fragmentation to transversely aligned quarkonium is slightly enhanced 
relative to longitudinally polarised. This net alignment corresponds to a ~ 5% asymmetry 
in the angular distribution of leptons produced in the subsequent decay of the quarkonium. 
We point out that the leading gluon contribution to tp and T production arises from 
the evolution of the charm quark fragmentation function, rather than from direct gluon 
fragmentation. 
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Quarkonium production in high energy processes is dominated by the fragmentation 
of heavy quarks and gluons For example, in Z° decay the short distance process 

Z° — > ifjgg is suppressed relative to the fragmentation process Z° — > i/jcc by a factor of 
m c/ m Z' corresponding to the small overlap of the if) wavefunction with a cc pair that has 
relative momentum of order mz- 

Recently it has been shown by Braaten and Yuan |3[] and Braaten, Cheung and Yuan 
that the process independent fragmentation functions for quarkonium production in 
high energy processes are calculable. For example, in Z° decay to a ijj of energy E 3> m^, 
the differential decay rate may be written in the factorised form 

dT(Z° -> ij;(E) + X) = J2 [ dzdt(Z° ->i(E/z)+X,(j) Di^{z,n). (1) 

Here dT(Z° — > i(E/z)+X, fi) is the differential decay rate for a Z° to decay to an on-shell 
parton of type i, while the fragmentation function D^^(z, fj.) gives the probability for the 
parton to split into a if) with momentum fraction z and another parton with momentum 
fraction 1 — z. The advantage of this factorised form is that the fragmentation function is 
independent of the production process. QCD corrections of the form log(mz /rn c ) neces- 
sitate the introduction of a factorisation scale /x; large logarithms of n/m c in Di^(z, fx) 
may be summed using the Altarelli-Parisi equations || @ . 

The fragmentation functions D c ^(z, (i) and Db^r(z,n) were calculated in 0], ex- 
plicitly summing over quarkonium polarisations. However, the fragmentation functions 
for transverse and longitudinally polarised quarkonium are individually calculable and are 
interesting in their own right. In this work we compute the fragmentation function for 
production of transversely polarised ift's from c quarks and T's from b quarks. 

The Feynman diagrams responsible for the fragmentation of a heavy quark into 
quarkonium are shown in fig. 1, where the black circle represents some quark-antiquark 
production process. For definiteness we consider Z° decay to i(jcc; however the fragmen- 
tation function we will derive is independent of the production process and is trivially 
extended to b fragmentation to T. Factorisation is only manifest in q ■ A = gauge (where 
q is the antiquark four- momentum) , in which the diagram in fig. 1(b) vanishes to leading 
order in m c /E. Following the manipulations in Eqs. (9)-(ll) of Ref. 0], we write the decay 
rate T(Z° %jj{e) cc), averaging over initial spins and summing over colours, as 

T(Z° - ^(e)cc) = J- / [dq] [dq] (2tt) - 1 - q) 

Z J (2) 

J 8nJ 2n \ z 1-zJ 1 K)l 
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where A(e) is the amplitude to produce a ip with polarisation e, pz is the Z° four- 
momentum, s = q 2 is the invariant mass of the original c quark, [dq] = d 3 q/(167r 3 qo), 
and ^ is the longitudinal momentum fraction of the ip. Choosing a frame such that 
q = (q°, 0, 0, g 3 ), we have z = (po + p^)/{qo + 93), where p is the four-momentum. 
Squaring the result from fig. 1 and using the identity 



to sum over transverse ip polarisations, we find 



Ei^)i 2 = jr^ DT ^ z) (-^ + ^)tr(r^), (4) 



where 



T . 256 2 2 J-2 
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2 y ~ \ z z 1 

^16ml(s - m 2 c )z 2 ^ 8(s - m 2 c ) 2 z 2 (l - z) 
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We have dropped terms suppressed by m c /E in our result, is the ifj decay constant, 
defined by 

(0\c^c\ij(p,e)) = Um^. (6) 

This is related to the nonrelativistic radial wavefunction at the origin, R^(0), by = 
Rip(0). Numerically, from the ip leptonic width, 410 MeV. The total decay 
rate to transversely aligned may now be written in the factorised form 



x fdz Hdsefs-^- -^U 1 d t ( 8 , z), 

Jo Jo V z 1 ~ Z J {s-m 2 c Y 



(7) 



where the decay rate to free quarks is 

r( Z° -> cc) = J- | [dg] [^](27r) 4 5( 4 )(p z - ? - g) 



X 



a /3 \ 

^ + PzPz t r(r Q |r^). 

m z j 
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Performing the integral over s in (0) , we find the expression for the fragmentation function 
to transversely aligned ifr's: 

D^(z,ix = 3m c ) =^^a s (3m c ) 2 /| 

(9) 



\~T^ — ~W ( 16 ~ 32z + 76z2 ~ 36zS + 6z *) ■ 



Adding the longitudinal polarisation, we recover the total fragmentation function D^^(z) 
given in ||], 

Dc+L( z ,n = 3m c ) = ^-ja s (3m c ) 2 tf 



2(1-2 



(10) 
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(16 - 322 + 722 2 - 322 3 + 52 4 ) . 



In experiments performed at a hadron collider, the gluon fragmentation function also 
plays an important role. The Altarelli-Parisi evolution of -D c _^ (z, fx) from the scale 
fj, = 3m c to a high energy scale fx — M (such as M = mz) will induce a gluon fragmentation 
function D g _^^ [z, M) via the gluon splitting function P g ^ C c(z). That is, a high energy 
gluon may split into a quark-antiquark pair of lower energy, either of which may then 
fragment to quarkonium. In the leading logarithmic approximation this effect, of order 
a.g ln(M/m c ), dominates over the direct contribution to gluon fragmentation computed in 
0, which is of order without the ln(M/m c ) enhancement. ^ 

It is equally straightforward to evolve the quark fragmentation functions (|^) and 
( |TUD to high energies using the Altarelli-Parisi equations and the quark splitting func- 
tions P c ^ cg (z) and P c ^ gc (z) = P c ^ cg (l — z). Because J Q dz P c -> cg (z) = 0, this evolu- 
tion softens the z distributions of D^^(z, fx) and D^+^z, jx) but leaves their integrals 

1 T ( T I T j *) 

Jo dz D c _X^ '(z,/x) unchanged ||. 

We now define £ to be the ratio of transverse to total fragmentation probabilities, 

As discussed in the previous paragraph, ( is independent of fx. Evaluating ( P"T| ) at fx = 3m c , 
we find ( = 0.69, to be compared with C = I f° r production of unaligned ^'s. Therefore 



t This applies only to ip and T production; gluon fragmentation to r/ c 's and ryt's is of order o? s , 
and so the direct gluon fragmentation calculated in pi dominates. 
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a small excess of transversely aligned ip's will be produced. Since the dependence on 
and m c drops out of ( |TT|) , the ratio £ is the same for ift's and T's. In leading logarithmic 
approximation, the corresponding ratio for gluon fragmentation to ift's has the same value 
and is \x independent, since the gluon fragmentation functions are induced only by quark 
fragmentation. Hence, at a hadron collider, where ip's are produced both by quark and 
gluon fragmentation, the fraction of transversely aligned ip's is also £ = 0.69. 

The asymmetry £ is measurable through the angular distribution of the leptons in the 
decay ijj — > £ + £~. Defining 9 to be the angle between the alignment axis and the lepton 
momentum, the angular distribution dT/dO in the if) rest frame has the form 

— K — — -ocl + acos 2 # (12) 

du 

where 

3^-2_ (13) 
2-C 1 ' 

£ = 0.69 corresponds to a = 0.053, about a 5% asymmetry. 

It is useful to compare this with the polarisation of ift's produced from nonleptonic 

b decay. Assuming the four-quark amplitude factorises (which should give a reasonable 

estimate, as the ip is a compact object on the hadronic scale Aqcd), a straightforward 

calculation of the decay b — > i/js — > £ + £~s gives 

2 2 

a = \ 1 w -0.46. (14 

This is in nice agreement with a recent CLEO inclusive measurement, which finds (see 
also g) 

a = -0.44 ±0.16 (15) 

(where we have added the systematic and statistical errors linearly). Since the inclusive 
branching ratio for b ijjX is (1.02 ± 0.05 ± 0.09)% 0, this decay may provide a sig- 
nificant background to if) production by fragmentation at high energy colliders. There is 
no comparable decay process to compete with b quark fragmentation for T production; 
therefore this mode should allow a cleaner measurement of £. 

We note that ip (or T) production can also arise from the fragmentation into P-wave 
quarkonium which subsequently decays electromagnetically to if) (or T), as mentioned in 
B. For weakly bound nonrelativistic quarkonium, the fragmentation into P-wave states 



is suppressed; however, in the charmonium system this may be an important additional 
source of ip's. 

In summary, we have calculated the leading contribution to the transverse fragmenta- 
tion function for quarkonium production from heavy quarks. We find that the and T are 
essentially unaligned, with only a ~ 5% angular anisotropy in the lepton distribution pro- 
duced in their electromagnetic decays. We also point out that the dominant contribution 
to gluon fragmentation into ip or T in high energy processes arises from Altarelli-Parisi 
evolution of heavy quark fragmentation functions, rather than from direct gluon fragmen- 
tation. 
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